AN EXPLICIT UNIVERSAL CYCLE FOR THE 
(n - l)-PERMUTATIONS OF AN n-SET 



FRANK RUSKEY AND AARON WILLIAMS 



Abstract. Wc show how to construct an explicit Hamilton cycle in the directed Cayley 
graph Cay({(T„, (T„_i} : Sji), where ct/c = (1 2 ••• k). The existence of such cycles was 
shown by Jackson (Discrete Mathematics, 149 (1996) 123-129) but the proof only shows 
that a certain directed graph is Eulerian, and Knuth (Volume 4 Fascicle 2, Generating All 
Tuples and Permutations (2005)) asks for an explicit construction. We show that a simple 
recursion describes our Hamilton cycle and that the cycle can be generated by an iterative 
algorithm that uses 0{n) space. Moreover, the algorithm produces each successive edge of 
the cycle in constant time; such algorithms are said to be loopless. 



1. Introduction and motivation 

There are many proofs in the mathematical hterature showing the existence of Hamilton 
cycles or Eulerian cycles in important families of graphs. However, turning these proofs into 
efficient algorithms often represents a significant challenge. 

An interesting case in point is the well-known De Bruijn cycle, which is a length k"' circular 
string over a /c-ary alphabet with the property that every length n string occurs as a sub- 
string. The existence of De Bruijn cycles is commonly presented in undergraduate discrete 
mathematics courses as a consequence of a certain graph being Eulerian. However, it is not 
widely known how to efficiently generate a De Bruijn cycle. In the authors' view two aspects 
of this question have particular importance. 

• Space, not time, is the primary enemy. A naive solution would be to build the 
graph and then use a Eulerian cycle algorithm to produce the cycle. This will be 
practical for small values of n and k, but for large values space will be the limiting 
factor long before time becomes a factor. In general, we need to be able to generate 
the Hamilton or Eulerian cycle without building the graph, or storing exponentially- 
long sublists. There are algorithms for building De Bruijn cycles that use space 0{n). 
The earliest of these is due to Fredricksen and Maiorana |3] and is presented in Knuth 

• The development of efficient algorithms reveals structure. It is often worth- 
while to turn a proof into an algorithm, or to develop an alternate proof, because the 
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process often results in a deeper structural understanding of the cycles being listed. 
For example, the efficient algorithm due to Fredricksen and Maiorana is based on 
necklaces, Lyndon words, and is related to pattern-matching and Lyndon factoriza- 
tion. 

As another example from the Hamiltom cycle domain. Fades, Hickey, and McKay [2] con- 
sidered the graph G{n, k) whose vertices are all length n bitstrings with density k and where 
two bitstrings are joined by an edge if they differ by transposing two adjacent bits. They 
showed that G{n, k) is Hamiltonian if and only if n is even and k is odd. The proof is 
inductive and relies on the fact that the graph has a spanning subgraph that is the prism 
of two "combs." However, it was not at all clear how to turn that proof into an efficient 
algorithm. Fventually an algorithm that mimics the proof was found that uses 0{n) space 
and take time 0(1) per bitstring generated [5]. 

In the present paper we are considering the construction of a "universal cycle" for the 
{n — l)-permutations of an n-set (which we take to be {1, 2, ... , n}). Here a universal cycle 
is a circular string of length n! what contains each of the n! different {n — l)-permutations 
as a (contiguous) substring. For example, 321312 is a such a universal cycle for n = 3, since 
its substrings are 32, 21, 13, 31, 12, and 23. 

More general universal cycles were introduced by Chung, Diaconis, and Graham [1] as a way 
of extending the de Bruijn cycle idea to combinatorial objects in general. The existence of a 
universal cycle for the fc-permutations of an n-set was shown by Jackson [6] when k < n. His 
proof sets up a certain natural Eulerian graph, call it Jk,n, and shows that any Fulerian cycle 
in that graph corresponds to the required universal cycle. However, no explicit construction 
of the cycle is indicated. The problem for k = n — 1 is discussed by Knuth [H] in Fxercise 
112 of Section 7.2.1.2. On page 121 of [8] we find the following quote: 

"At least one of these cycles must almost surely be easy to describe and to 
compute, as we did for de Bruijn cycles in Section 7.2.1.1. But no simple 
construction has yet been found." 

The purpose of this paper is to provide such a description and computational method. We 
will show how to construct a particular universal cycle. Our algorithm takes space 0{n) and 
uses a constant amount of time between successive outputs of characters in the cycle. To be 
precise regrading the space requirement: The algorithm uses a constant number of arrays, 
each with 0{n) indices, and each storing integers of value at most n. Similarly regarding 
time, we use a constant number of operations (comparisons, increments, decrements, and 
parity tests) on integers of value at most n. 

Universal cycles for the permutations of an n-set are not directly possible unless n < 2. 
However, every (n — l)-permutation of an n-set can be uniquely extended to a permutation 
of an n-set by appending the unique missing symbol. Thus, universal cycles for (n — 1)- 
permutations can be viewed as universal cycles for permutations. For example, 321312 
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produces the permutations 321, 213, 132, 312, 123, and 231, where the appended missing 
symbols are underhned. For this reason, our results add to the already sizeable literature 
on generating permutations. A good survey is provided by Sedgewick [11] and more recent 
developments are to be found in Knuth [8]. 

We don't expect our algorithms to be a fast way to generate permutations using the usual 
model of computation, since at least n — 1 of the n values change at each step. However, 
they will be fast if a circular representation is used; for example, when using linked lists or a 
circular array. In a circular array we maintain a start position and do arithmetic on indices 
mod n. They will also be fast if the permutation is stored as a computer word. For example, 
we can store the permutations up to n = 16 by dividing 64 bit words into 4 half-bytes. The 
shifts can then be accommodated in a few machine instructions. 

Finally, we mention that additional symbols can also be used to create universal cycles 
whose substrings are order isomorphic to permutations. For example, 421423 produces the 
permutations 321, 213, 132, 312, 123, and 231. Recently Johnson [7] proved a conjecture 
in [1] by showing that n + 1 symbols are always sufficient for constructing these universal 
cycles. 

The paper is organized as follows. In Section 2 we give our explicit construction as a certain 
recursively defined string. Then, in Section 3, we show that this string can be generated 
by an algorithm that uses only a constant amount of computation between the output of 
successive symbols of the string — the first such algorithm for a universal cycle. In Section 
4, we give further properties of our recursive construction; first some results on the number 
of (Tn or cTn_i Operations that are used, then that our ordering has an efficiently computable 
ranking function, and finally that it is "multiversal," in a sense to be described later. We 
conclude with Section 5, which contains some open problems. 

2. An explicit construction 

Initially, we will couch our discussion in terms of finding Hamilton paths in certain di- 

— > 

rected Cayley graphs. Cayley graphs are denoted X = Cay({ai, a2, ■ ■ ■ , «fe} : G). Here 
{ai,a2, ■ ■ ■ ,afc} is a generating set of a group G. The vertices of X are the elements of 
G and the edges are all of the form g aig; these edges are usually thought of as being 
labelled with In an undirected Cayley graph, if a is in the generating set, then its inverse 
a~ is also in the generating set. Driven by the question of Lovasz of whether there is a 
Hamilton cycle in all undirected Cayley graphs, there is a significant literature of results 
about Hamilton cycles in Cayley graphs. A survey may be found in Gallian and Witte [1]; 
see also Pak and Rados Radoicic |10j . 

In the solution to Exercise 112 of Section 7.2.1.2 Don Knuth implicitly poses the problem 
of finding an explicit expression for universal cycles of [n — l)-permutations of an n-set [8]. 
This problem is equivalent to generating permutations of an n-set by rotations of the form 
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(12 • ■ ■ n) or (1 2 ■ • • n— 1); i.e., it is equivalent to asking whether the Cayley graph 
— > 

S„ := Cay({(j„, cr„_i} : §„) 

is Hamiltonian. We use E>k,n to denote the set of /c-permutations of the n-set [n] = {1,2, ... ,n}. 
In the case where k = n we use S„. Although we do not use this fact below, it is interesting 
to note that a short proof reveals that the graph S„ is the line graph of the Jackson graph 

J n— l,n- 

Consider the binary string Sn defined by the following recursive rules. The base case is 
5*2 = 00. Let Sn = X1X2 ■ ■ ■ Xn\ where x denotes flipping the bit x. Then, for n > 2, 

(1) Sn+i := 001"-^ Xi001"-2 X2 ■ ■ ■ 001"-2 Xn\. 

We use above the usual convention that if w; is a string and m is an integer then w"^ is w 
concatenated together m times, w"^ = ww ■ ■ -w; also is the empty string. 

Below we list S3, 6*4, and 6*5. Each Si is of the form ww since 00 has this property and the 
recurrence ([1]) preserves it. 

^3 = 00 00 = 00 1 00 1. 

^4 = 001 001 001 1 001 001 001 1 = OOl l OOl l OOl O OOl 1 001 1 001 0. 
^5 = ( 0011 1 0011 1 0011 0011 0011 1 0011 1 

0011 0011 0011 1 0011 1 0011 0011 1 f 



Now define the mapping by — (J„, and 1 where 0"^ = (1 2 ■ ■ ■ k). 

Theorem 2.1. The list (p{Sn) is a Hamilton cycle in the directed Cayley graph 

Proof. In listing the Hamilton cycle we use one-line notation for the permutations, starting 
with n n—1 ■ ■ ■ 2 1, and think of the cycles cr„_i and an as acting on the positions in the 
one-line notation. Thus, in a slight abuse of notation, 

(piSs) = 321, 213, 132, 312, 123, 231, 

since S3 implies the successive application of a^, a^, 02, 03, 03, and finally Oi to map the 
last permutation to the first. 

Our proof strategy is to give an explicit listing of permutations of [n] with the required 
properties and then show that it is equivalent to ([1]). Recursively define a circular list 
n(n) = 11(72)0, n(r;,)i, n(?7.)„!_i of permutations of [n]. For small values of n, define 
n(l) = 1, n(2) = 21, 12, and 11(3) = ^(S's). Every n-th permutation of n(n) is defined as 
follows. 



(2) n(n)j>, := nn(n - 1)^-. 
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The n — 1 permutations that follow mr, where n = Il{n — are defined to be 

(3) (Tnimr), alinn), a„_i(o-^(n7r)), <Zi(a2(n7r)). 

The list n(4) is shown in Tabled column (d). The permutation nvr followed by the permuta- 
tions above comprise the sublist n(n)j„, n(n)j„+i, . . . , n(n)(j+i)„_i and these permutations 
are all distinct since the position of n is successively in the n different positions 1, n, n — 1, 
. . . , 2. Furthermore, because we can recover tt from any permutation in this sublist, the 
uniqueness of every permutation in Il{n) follows inductively from the uniqueness of every 
permutation in n(n — 1). 

It remains only to prove that successive permutations differ by cr„ or cr„_i and that the list 
is circular. It is clear from ([3]) that successive permutations differ by (T„ or an-i, except 
for those that precede the one of the form nn successively followed by nir. Let arz be a 
permutation of 1,2, ... ,n — 1 where a and z are numbers and r is a sequence (of length 
n — 3). Note that the last permutation of ([3]) is 

(^n-l{(^l{narz)) = a;^\{al{naTz)) = a~\{Tzna) = znra. 

Now suppose that vr = n(n— 1)^ = arz and vr' = n(n— Inductively, either vr' = cr„_i(7r) 
or vr' = o"„_2(7r). Observe that 

(4) (Tn{znTa) = nraz = n(T„_2(ar2;), and 

(5) an~i{znTa) = nrza = nan-i{aTz). 

Since the successor of arz is either cr„_2(«T-2) or cr„_i(ar2), the transition to the permutation 
vr' is also of the correct form; successive permutations in n(n) differ by cr„ or cr„_i. The 
circularity of the list follows inductively from the circularity of the list n(n — 1) (alternatively 
we could use Lemma 12.21 below). Furthermore, in terms of the mapping defined earlier, 
the bits are fiipped; a (cr„_i) transition in n(n — 1) becomes a 1 (o"„_i) transition in n(?7,) 
by ([5]), and a 1 ((7„_2) transition in liiji — 1) becomes a ((T„) transition in n(n) by (jl]). □ 

Lemma 2.2. Any Hamilton path in S„ is, in fact, a Hamilton cycle. 

Proof. Suppose that 11 = Hi, 112, ... , H^i is a Hamilton path in S„ that is not a Hamilton 
cycle. In particular (T„(H„,!) ^ Hi and cr„,_i(H„!) ^ Hi. Thus H„! ^ a" (Hi) and H„! ^ 
(j~_i(Hi). We must then have that cr~(Hi) — > (T„_i((T~(Hi)) and (j~_i(Hi) — > cr„(cr~_i(Hi)) 
are distinct edges in n(n). However, an easy calculation shows that a^^^an = cr^an-i = 
(n— 1 n) and thus these permutations are identical. This contradiction shows that H is a 
Hamilton cycle. □ 

The proof shows that in fact the lemma is true for any Cayley graph on two generators p 
and r for which T~p is an involution. 

The universal cycle for (n — l)-permutations of [n] is obtained by recording the first symbol 
in each of the permutations in H(n). We use Un to denote the resulting universal cycle. 
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3. A LOOPFREE ALGORITHM 

Suppose that in our recurrence ([1]) for 5'„+i that for each "new" bit we record the value 
n, and apply this idea recursively. Call the corresponding new sequence Rn+i- That is, 
i?2 = 11; and for n > 1, 

Rn+i = n yin y2---n y^], 

where Rn = yiy2 ■ ■ ■ Vni- For example 

i?4 = 333 2 333 2 333 1 333 2 333 2 333 1. 

The sequence R4 is exactly the sequence that is obtained by recording the most significant 
position that changes when counting with the multi-radix numbers with parameters 2x3x4, 
when the numbers are indexed 1, 2, 3, from left-to- right. See Tabled! columns (a) and (b). 
In general, Rn gives us the positions when counting with multi-radix numbers 2 x 3 x • ■ ■ x n. 

These observations suggest that we may be able to efficiently generate the Sn sequence 
by modifying the classic algorithm for counting with multi-radix numbers. In the classic 
algorithm the multi-radix number is stored in the array a„_i ■ ■ • 0201 and j is used to represent 
the rightmost, or smallest, index where aj is not at its maximum value. The next multi- 
radix number is obtained by incrementing aj and setting all values to its right to 0. Now 
suppose that we just incremented the integer in position j so that the multi-radix number 
IS fln— 1 ■ ■ ■ CLjClj — i ■ ■ ■ Oil = CLn—t ' ' ' CL ,0 • • • 0. Then the corresponding R^ value is n — j and so 
the non-recursive part of the 5*^ sequence that we are listing is going through the pattern 
QQX'^-J-i or the pattern 110""-'"^, depending on whether j is odd or even, respectively. For 
proposition P we use the notation \P\ to mean the value 1 if P is true and the value if 
P is false; also © denotes exclusive-or. The expression \j even © aj < 1] gives the correct 
value of the bit to be output. Below is the entire algorithm, rendered in pseudo-code. 



a„+ia„, ■ ■ ■ ai ^ ■ ■ ■ 


0; 


repeat 








while aj = n — 


j do aj ^0; j ^ j + 1; od; 

< 11 ); 


output( [j even 


aj <— aj + 1; 




until j > n; 





There is an loopless algorithm for listing multi-radix numbers as a Gray code in which the 
value in only one position changes and that change is by ±1 (see, for example, Williamson 
[12], pg. 112, or Knuth [8], pg. 20). Together with the ideas used in the previous "counting" 
algorithm, we can adapt those loopless algorithms to get a loopless algorithm for generating 
Sn or our universal cycle. In the Gray code for multi-radix numbers, the values in a given 
position alternately increase and decrease. Furthermore, the values change in exactly the 
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(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


(g) 


234 




234 


54 


§4 


Ua 


rank 


000 


3 


000 


. . 


4321 


4 





001 


3 


001 


. . 


3214 


3 


1 


002 


3 


002 


. . 1 


2143 


2 


2 


003 


2 


003 


. 1 . 


1423 


1 


3 


010 


3 


013 


. . 


4213 


4 


4 


Oil 


3 


012 


. . 


2134 


2 


5 


012 


3 


Oil 


. . 1 


1342 


1 


6 


013 


2 


010 


. 1 . 


3412 


3 


7 


020 


3 


020 


. . 


4132 


4 


8 


021 


3 


021 


. . 


1324 


1 


9 


022 


3 


022 


. . 1 


3241 


3 


10 


023 


1 


023 


. . 


2431 


2 


11 


100 


3 


123 


. . 


4312 


4 


12 


101 


3 


122 


. . 


3124 


3 


13 


102 


3 


121 


. . 1 


1243 


1 


14 


103 


2 


120 


1 


2413 


2 




110 


3 


110 


. . 


4123 


4 


16 


111 


3 


111 


. . 


1234 


1 


17 


112 


3 


112 


. . 1 


2341 


2 


18 


113 


2 


113 


. 1 . 


3421 


3 


19 


120 


3 


103 


. . 


4231 


4 


20 


121 


3 


102 


. . 


2314 


2 


21 


122 


3 


101 


. . 1 


3142 


3 


22 


123 


1 


100 


. . 


1432 


1 


23 



Table 1. (a) Counting in multi-radix base 2x3x4, (b) the sequence, (c) 
the corresponding multi-radix Gray code, (d) indented version of 5*4, (e) the 
list n(4), (f) the universal cycle U^, and (g) the rank of each permutation. 



positions given by the Rn sequence. In the implementation we maintain a direction array 
d where -|-1 means increase and —1 means decrease We also maintain an array / of "focus 
pointers" which allow instant access to the next position whose value will change (we set 
fn = n+1 (instead of n) so that the last iteration is handled correctly). See Table [1], column 
(c), for an example. 

Thus the values of j from the counting algorithm are exactly the same in the Gray code 
algorithm, except that in the Gray code algorithm j is the position where a value changes. 
The only complication arises because the values in a given position can be decreasing, and 
so the test "a^ < 1" is not sufficient. Fortunately, all algorithms that looplessly implement 
the Gray code maintain an array of directions dn-i ■ ■ ■ d2di for each position, where di G 
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{+1, —1}, indicating whether the values in that position are currently increasing (+1) or 
decreasing (—1). li dj = +1 then we can continue to test aj < 1, but to account for dj = — 1, 
we need to test 

{ttj < 1 and dj = 1) or (a^ > n — j — 1 and dj = — 1). 

We can "optimize" this condition. Notice that the test (oj < 1 and dj = 1) can be replaced 
by aj — dj < 0. This change is possible because if dj = —1 then aj — dj is guaranteed to be 
greater than zero because aj > 0. Therefore, if aj — dj < 0, then this immediately implies 
that dj = 1 and so % — 1 < 0, which is equivalent to the original test aj < 1. Likewise, the 
test {aj > n — j — 1 and dj = —1) can be replaced by aj — dj > n—j. Below is our loopless 
algorithm in full detail. 



Q'n+lO'n ■ 


■■ai ^ ■ ■ ■ 


0; 




d-adn—l 


■ ■ rfi ^ 1 1 1 ■ ■ ■ 


1; 




fnfn—1 


■ ■ /i ^ n+1 n—1 


n-2 ■•■1; 




repeat 








J 


^ /i; fi - 1; 






output( 1 j even © 


[aj -dj<0 


or aj - dj >n-j)j); 


aj 


^ aj + dj] 






if 


aj = or aj = n 


— j then dj 


' — df, fj ^ fj+i, f]+i ^ J + 1; fi; 


until j 


> n; 







It is also possible to output the universal cycle itself in a loopless manner, but an additional 
circular array is required to hold the current permutation. To follow are the details. Define 
an array 7ri7r2 ■ ■ ■ vr„ initialized to n n—1 ■ ■ ■ 1 and an index t that will be incremented mod 
n on each iteration of the algorithm. We will think of vr as a circular array. The index t is 
the position of the last element of tt, so initially t = n. As each bit of Sn is determined, we 
will ouput the first element of tt (i.e., the one in position t + 1). If the bit is a 1, so that 
(Tn-i is acting on vr then we need to swap the last two elements: VTt-i ^ irt. In other words 
the output statements in the preceding code fragments is replaced with the following code 
where expr is the expression inside of the output statement in either the previous counting 
algorithm of the previous loopless algorithm. 



t' ^t] t ^ (t+l) mod n; 

OUtput( 71(1+1) modn ); 

if expr = 1 then ttj ^ tt^/ fi; 



Finally, we note that every permutation can be output in a circular fashion by outputting vr 
and t. We could also use a linked list, which would give a loopless permutation generation 
algorithm. 
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4. Further properties 



In this section we explore further properties of S„ and our Hamilton cycle. 

4.1. How many of each rotation is used? It is clear from the recurrence relation ([1]) 
that the number, call it /„, of (T.„'s in (p{Sn) satisfies the recurrence relation 



Interestingly, this sequence appears in OEIS [S] as Al22972(n + 1) as the solution to the 
"symmetric" recurrence relation a{n + 1) = (n — 1) ■ a{n) + n ■ a{n — 1). The values of /„ for 
n = 1..10 are 1, 2, 4, 14, 58, 302, 1858, 13262, 107698, 980942. 

Consider the cosets induced by an', there are n\/n = (n — 1)! of them. In a Hamilton 
cycle there must be at least one an-i edge that leaves each coset, and thus there must be 
at least (n — 1)! of them. Alternatively, consider the cosets induced by an-i', there are 
n!/(n — 1) = n ■ {n — 2)! of them. In a Hamilton cycle there must be at least one (t„ edge 
that leaves the coset, and thus there must be at least n ■ (n — 2)! of them. We can make a 
stronger statement regarding the (T„ edges. 

Lemma 4.1. The least number of an edges in any Hamilton cycle in S„ is 2n{n — 2)! — 2. 

Proof. First, observe that 

aZiana^_ian = (n— 1 n){n—l n) = id. 

The two cr„ edges above are incident with the same unordered pair of cosets induced by 
an-i- Thus if we contract each coset into a singe super- vertex, then the resulting graph, 
call it Qn, is undirected in the sense that every directed edge is paired with an edge in a 
2-cycle. Furthermore, it is not hard to see that if one of those (j„ edges is used in a Hamilton 
cycle, then so must the other. Thus a Hamilton cycle in X„ becomes a connected spanning 
subgraph of Qn- Since a minimal connected spanning subgraph is a spanning tree, and any 
spanning tree has n ■ (n — 2)! — 1 edges, the number of o"„ edges is at least 2n{n — 2)! — 2. □ 

Figure [1] shows the Cayley graph X^. Note that the contracted graph is the 3-cube. The 
red edges show the Hamilton cycle 5*4. In this case Sn corresponds to a spanning tree in Qn, 
but this is not the case for n > 6. 




This recurrence relation can be iterated to obtain 



n-l 



/„ = 2(-ir-3^(-l)'=(n-A;)!, 



k=l 

from which it follows that 



/„, ~ 3{n - 1)! or ^ ~ -. 

nl n 
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V 



123 




234 





Figure 1. The Cayley graph S4 on the left. The graph Q4 on the right. The 
thick (red) edges indicate the Hamilton cycle S4. 



4.2. Ranking. The rank of a permutation vr is the value r for which Il(n)r = vr. Our 
recursive equation for the rank depends on the position of n within the permutation being 
ranked. From the definition of Il(n) we can infer that 



rank(aia2 ■ ■ ■ ttk-iuak+i ■ ■ ■ an) 

'0 ifn = l, 

n ■ rank(a2a3 ■ ■ ■ a„) if /c = 1, 

n — k + 1 + n ■ rank(a„afc+i ■ ■ ■ a„_iai ■ ■ ■ a^) if A; > 1. 

The expression n — k + 1 accounts for the position of the n, and the rest comes from the 
recursive part of the definition of n(?T,). We can also express the rank as 



{0 if a = /3 = e, 

n ■ rank( /5 ) if a = e, 

n — \a\ + n ■ rank( a{P)a ) otherwise , 

where cr(/?) is f3 rotated one position to the right. 

Implemented in the obvious manner, these recurrence relations lead to algorithms that use 
O(n^) arithmetic operations on integers as large as n\. 

4.3. Multiversal Cycle Property. In this section we prove that the sequence n(?T,) = 
n(n)o, n(?T,)i, . . . , n(n)„!, written out as a long string of symbols by concatenating each 
permutation, is a "multiversal cycle". We denote this "flattening" of n(n) as For 
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example, consider ]J(3) = 321 213 132 312 123 231. Starting in positions 0,1, or 2 and 
advancing the position in increments of 3, recording the first two symbols, we obtain 






32 


21 


13 


31 


12 


23 


1 


21 


13 


32 


12 


23 


31 


2 


12 


31 


23 


21 


32 


13 



In each case a complete set of all 2-permutations of [3] is obtained. The purpose of this 
section is to prove that this property holds in general. 

Definition 4.2. A multiversal cycle for the (n — l)-permutations of an n-set is a circular 
string a^ai ■ ■ ■ a^^i of length N = n ■ n\ such that, for all m = 0, 1, . . . , n — 1, 

(7) \^Clm+in ' ' ' Oim+in+n—2 | ^ 1)2,..., 77.!} l,n) 

where arithmetic in the indices is taken mod n. 

Before getting to the main theorem in this section we prove a technical lemma. 
Lemma 4.3. For all i 7^ 0, — 1 mod n, if = a^ai ■ ■ ■ a^-i, then 

(^i (^i-\-n—l- 



Proof. Because i ^ mod n the numbers and aj+„_i lie in two successive permutations of 
n(n). The conclusion now follows since successive permutations differ by (T„ or cr„_i. The 
i ^ —1 mod n condition is necessary when they differ by Cn-i- □ 

Theorem 4.4. The string = aoOi ' ' ' Oni-i ^-5 o- multiversal cycle. 

Proof. The proof is by induction on the value m in the definition. The base case m = sat- 
isfies because n(n)o, n(n)i, . . . , n(n)n!-i is a listing of all permutations of [n], so ignoring 
the last character of each permutation gives a complete listing of all (n — 1) -permutations 
of [n]. Similarly, when m = 1, ignoring the first character of each permutation also gives a 
complete listing of all (n — 1) -permutations of [n]. We now argue by contradiction. Suppose 
that there are some values m > 1, i and i', with i 7^ i', such that 

(8) OYn+in ' ' ' Oim+in+n~-2 ^m+i'n ' ' ' (^m+i'n+n—2- 

Inductively, we know that 

(^m—l+in ' ' ' (^m~l+in+n~2 7^ '^m— 1+j'n ' ' ' (^m—l+i'n+n—2- 

Thus it must be the case that am-i+m 7^ dm-i+i'n- However, applying Lemma 1473) to am-i+m 
and am-i+i'n gives 

(^m—l+in (^m—l+in+n—1 Sud dm— 1+i'n 0,m—l+i'n+n—l} 
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SO long as m 7^ 0, 1. But by ([H]) we now have 

(^m—l+in '^m— 1+m+n— 1 (^•m+in+n—2 '^m+i'n+n— 2 (^m—l+i'n+n—l '^m— l+i'n; 

which is a contradiction. □ 

The careful reader will have noted that Lemma 14.31 and Theorem 14.41 apply to any Hamilton 
cycle in S„ since the only property that we use is that successive permutations differ by (j„„i 



5. Final Remarks, Open Problems 

In this paper we have developed an explicit algorithm for generating a universal cycle for 
the (n — l)-permutations of an ra-set. This is the first universal cycle for which a loopless 
algorithm has been discovered. 

Below is a list of open problems inspired by this work. 

• Can the results of this paper be extended to /c-permutations of [n] for 1 < k < n — 17 

• Among all Hamilton cycles in S„ we determined in Lemma 14.11 the least number of 
(T„ edges that need to be used in a Hamilton cycle in S„. What is the least number 
of (T„_i edges that need be used? In our construction, the number of (T„ edges is 
asymptotic to 3/n and the number of (T„_i edges is asymptotic to {n — 3)/n. Is there 
a general construction that uses more (T„ edges than edges? 

• Can the results of this paper be extended to the permutations of a multiset? That 
is, given multiplicities no, ni, . . . ,nt, where rij is the number of times i occurs in the 
multiset and n = tiq + ni + ■ ■ ■ + rit, is there a circular string Cli(l2 ■ ■ ■ CLn 

of length 

= ( ^ ) with the property that 

\no,ni,...,nt/ f f j 

{tti Cj+i ■ ■ ■ ai+n-2 i(oj, Oj+1, • • • , ai+n-2) I'i- < i < N} 

is equal to the set of all permutations of the multiset. Since the length of ajaj+i ■ ■ ■ aj+„_2 
is n — 1 it is not a permutation of the multiset; one character is missing. The function 
L gives the missing character. We call these strings shorthand universal cycles. The 
current paper gave a shorthand cycle for permutations of [n]. 

• It would be interesting to gain more insight in to the ranking process. Is there a way 
to iterate the recursion so that it can be expressed as a sum? 
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